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Abstract 

We derive the generalization of Wigner's causality bounds and Bethe's integral formula for 
the effective range parameter to arbitrary dimension and arbitrary angular momentum. We also 
discuss the impact of these constraints on the separation of low- and high-momentum scales and 
universality in low-energy scattering. Some of our results were summarized earlier in a letter 
publication. In this work, we present full derivations and several detailed examples. 
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I. INTRODUCTION 



Causality plays a fundamental role in physics. The principle that no action can be ob- 
served before its cause puts important constraints on physical theories. In classical electrody- 
namics, causality leads to the Kramers-Kronig relations which relate the real and imaginary 
parts of the dielectric constant. In quantum mechanics, causality requires that no scattered 
wave is produced before the incident wave first reaches the scatterer. For finite-range inter- 
actions the constraints of causality on elastic phase shifts were first investigated by Wigner 
[1]. To illustrate the underlying physics, we consider a wavepacket of outgoing spherical 
waves in d spatial dimensions, 

POO _ 

/out(r)= / rfpe^P7out(p), (1) 

where normalization factors and the r"*-'^"^-'/^ radial dependence are absorbed into the def- 
inition of /out('")- When this wavepacket is scattered, the S-matrix multiplies asymptotic 
outgoing states by a phase factor e^*^*^^\ where 5{p) is the elastic phase shift. We assume the 
momentum distribution /out(p) is sharply peaked around some nonzero value p. If /out(^) is 
the scattered wavepacket, then 



Jo 



As a consequence, the wavepacket is shifted forward by Ar = —26'{p) relative to the 
wavepacket with no scattering. If we consider the wavepacket as a function of time, the 
same shift can be interpreted as a time shift or delay for the scattered wavepacket [1], 



d6 

At = 2 — 
dE 



(3) 



where E is the energy corresponding with p. The radius shift and time delay are sketched 
in Fig. [Tj A classical analysis of particle trajectories suggests that if the interactions have 
a finite range R, then causality requires —S'{p) < R. While this argument is qualitatively 
correct, it ignores the quantum mechanical spread of the wavepacket in space. The precise 
causality bound and the consequences of the bound are the subject of the present analysis. 

In a previous letter, we have explored the impact of causality constraints on low-energy 
universality [2j . Low-energy universality can appear when there is a large separation between 
the short-distance scale of the interaction and the long-distance scales relevant to the physical 
system. One example of low-energy universality is the unitarity limit. The strict unitarity 
limit refers to an idealized system where the range of the interaction is zero and the S'-wave 
scattering length is infinite. In experiment, the strict unitarity limit can not be reached 
because real systems always have finite range interactions. Since finite range effects are 
suppressed at low energies, one also refers to systems in the unitarity limit if they have 
infinite scattering length but finite range. In nuclear physics, cold dilute neutron matter 
is close to the unitarity limit since the neutron-neutron scattering length is much larger 
than the range. Stable many-body systems with infinite scattering length can be created 
in experiments with two types of ultracold fermionic atoms using Feshbach resonances. 
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FIG. 1: Time versus radial distance. The scattered wavepacket is shifted in distance by —2S'{p) 
and shifted in time by 2 J^j^.. 

The implications of causality for universality in 5'-wave scattering are well understood 
P]. In this case, the unitarity limit can always be reached by tuning the scattering length 
to infinity. For reviews of recent cold atom experiments exploring physics of the unitarity 
limit, see Refs. |ll[5]. Theoretical overviews of ultracold Fermi gases close to the unitarity 
limit and their numerical simulations are given in |6l [7] . A general review of universality at 
large scattering length can be found in [8]. 

The implications of causality for large scattering length physics in higher partial waves are 
more intricate. Because of causality, it is not always possible to reach the unitarity limit by 
tuning external parameters [21 19] . Several experiments have investigated strongly-interacting 
P-wave Feshbach resonances in ^Li and "^^K pUl - fn] . A key question is whether the physics 
of these strongly-interacting P-wave systems is universal, and if so, what are the relevant 
low-energy parameters. A positive answer to this question would provide a connection 
between the atomic physics of P-wave Feshbach resonances and the nuclear physics of P- 
wave alpha-neutron interactions in halo nuclei. Some progress in addressing these questions 
has been made with low-energy models of P-wave atomic interactions [T51420] and P-wave 
alpha-neutron interactions [2T] - [2^ . A renormalization group study showed that scattering 
should be weak in higher partial waves unless there is a fine tuning of multiple parameters 
[9]. Complementary work was carried out by Ruiz Arriola and collaborators. A discussion 
of the Wigner bound in the context of chiral two-pion exchange can be found in [23] while 
correlations between the scattering length and effective range related to the Wigner bound 
were discussed in [2B] . 

In Ref . [2] , we have addressed the question of universality and provided expressions for the 
causality constraints in arbitrary dimension d and for arbitrary angular momentum L. Our 
analysis applies to any finite-range interaction that is energy independent, non-singular, and 
spin independent. Our results can be viewed as a generalization of the analysis of Phillips 
and Cohen [3], who derived a Wigner bound for the S'-wave effective range for short-range 
interactions in three dimensions. In the current paper, we present full derivations of Wigner's 
causality bounds fT] and Bethe's integral formula [21] for the effective range parameter to 
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arbitrary dimension d and angular momentum L. The extension of Bethe's integral formula 
for d = 3 and L > was first derived by Madsen [28]. We discuss the impact of these 
constraints on the separation of low- and high-momentum scales and universality in low- 
energy scattering using several detailed examples. 

The paper is organized as follows. In Sec. [Tl| we set the stage by defining angular mo- 
mentum, scattering phase shifts, and the effective range expansion for the general case of 
d spatial dimensions. The equation for the radial wave function and the Wronskian iden- 



tity for two solutions with different energies are derived in Sec. Ill This identity is used in 



Sec. |lV]to derive the causality bound on the effective range for zero energy and the general 
bound for finite energies. In Sec. [V], we discuss the impact of these bounds on low-energy 
universality. In particular, we address the question of the unitarity limit. The impact of the 
causality bounds is illustrated in detail in Sec. VI using three different examples: a spherical 
step potential in d dimensions, alpha-neutron scattering which corresponds to exponentially 
bounded interactions, and the long-range van der Waals interaction. We end with a sum- 
mary and outlook in Sec. VII[ Finally, in the Appendices we explicitly check the Wronskian 
identities for all possible values of d and L and demonstrate the equivalence of our causality 
bound with Wigner's original bound on the energy derivative of the phase shift. 



II. PRELIMINARIES 

Our goal is to generalize Bethe's integral formula for the effective range parameter and 
the related causality bound for arbitrary d and L. We start with some general definitions. 

We consider two non-relativistic spinless particles in d dimensions with a rotationally- 
invariant two-body interaction. The interaction is assumed to be energy independent and 
have a finite range R beyond which the particles are non-interacting. Let fi be the reduced 
mass and p^/ (2/i) be the total energy. For d > 1 angular momentum is specified hj d — 1 
integer quantum numbers [211 EO] 

L={Mi,--- ,M,_i}, (4) 

satisfying 

|Mi I < M2 < ■ ■ ■ < Md-2 < Md-i. (5) 

We let L label the absolute value of the top-level quantum number, |M(^_i|. For example 
when c? = 3, Mi is the angular momentum projection M = —L, . . . ,L and M2 = L = 
0, 1, 2, . . . is the total angular momentum. In = 2, there is only one angular momentum 
quantum number \Mi\ = L = 0, 1, 2, . . .. The case of one spatial dimension is special since 
continuous rotations do not exist. Instead of rotational invariance, the key symmetry in 
one dimension is invariance under parity. We assume a parity-symmetric interaction and 
write L =L = for even parity and L =L = 1 for odd parity. In the following all results 
we derive for rotationally-invariant interactions in d > 1 are also valid for parity-symmetric 
interactions in d = 1. 

We analyze the two-body system in the center-of-mass frame using units with h = 1 
for convenience. The full wave function for reduced mass n and energy p^/(2/i), can be 
separated into a radial part and an angular part via 

^S(r) = RiUr)Y^n (6) 
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where the Yi^l-r) are hyperspherical harmonics. The hyperspherical harmonics in d dimen- 
sions satisfy the orthogonahty condition, 

J Y*{f)Yi^,{f) dnd = 6-L,,i., (7) 

where is the sohd angle and the sum rule [29] . 

2^ Yi^{r)YL{r) = C/ (r ■ r ). (8) 



L, L fixed 



Here ^"^ is a Gegenbauer polynomial, and the normalization factor A^^ is given by, 

(c/-2)!!27r^/2 

Nd = 77-, • (9) 

For d = 2 the expressions are defined in the limit — )■ 2. 

We can check that these give the expected sum rules for d < 3. In one dimension we get 
even and odd parity functions of f ■ f', 

{d + 2L- 2) {d - A)\\cT-\r ■ r') = { ^ \^^l^L = l. ^^^^ 



For two dimensions, we have a sum of modes e^*^^ for 9 = cos " ir ■ r' 



lim(d + 21 - 2) (d- 4)ncf -(f . f) = I |i J-l7,^,0,°fo, i > 0. (") 
In three dimensions, we recover the Legendre polynomials, 

{d + 2L- 2) {d - 4)!!^^ 2"^(f ■ f') = {2L + l)Pi(f ■ f'). (12) 

The scattering phase shifts are directly related to the elastic scattering amplitude fL,d{p), 
where 

p2L 

In addition to having finite range, we assume also that p'^^~^'^~'^ cot 6 L,d{p) does not diverge 
at p = and the interaction is not too singular at short distances. Specifically, we require 



that the effective range expansion defined below in Eq. ( 15 ) converges for sufficiently small 
p. Moreover, the reduced radial wave function m^^^ defined in Eq. (21) must satisfy that 

^u^l\ is finite and u^^\ vanishes as r — )■ 0. In [30], these short-distance regularity conditions 
are shown to be fulfilled for interactions arising from a static potential, 

W{r,r') = V{r)6{r -r'), (14) 

provided that V{r) = 0(r~^+'^) as r — for positive e. In our discussion, however, we make 
no assumption that the interactions correspond to a local potential. 
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The effective range expansion for general dimension d and angular momentum L is 



2L+d-2 



cot 6L,d{p) - 5(dmod2),0- In {pPL,d) 



TC 



(15) 



n=0 



The term 5(dmod2),o is for odd d and 1 for even d. ai^d is the scattering parameter/ rL^d 
is the effective range parameter, and V^l are the n*^-order shape parameters, pi^ is an 
arbitrary length scale that can be scaled to any nonzero value. The rescaling results in a 
shift of the dimensionless coefficient of p2L+d-2 ^^^q right-hand side of Eq. (15), and we 
define pi^d as the special value for pL^d where this coefficient is zero. 

Throughout our discussion we use the examples of S-wave and P-wave scattering in three 
dimensions to illustrate general formulas. For d = 3 and L = the scattering amplitude is 



/o,3(p) oc 



pcot5o,3(p) - ip' 



and the effective range expansion is 

pcot5o,3(p) = -— + Jro,3P^ + E(-1)"^'^&V"^' 
For d = 3 and L = 1 we have 



0-0,3 



n=0 



/1,3(P) oc 



and 



cot 5i,3(p) 



ai,3 



p^ cot 5i,3(p) — ip^ 

oo 



(16) 



(17) 



(18) 



(19) 



n=0 



III. RADIAL EQUATION AND WRONSKIAN IDENTITY 

The next step is the derivation of the Wronskian identity for the solutions of the radial 
Schrodinger equation. The causality bound then follows directly from this identity. 

The interaction is assumed to have finite range R beyond which the particles are non- 
interacting. With the interaction written as a real symmetric operator with kernel W{r, r'), 
the radial Schrodinger equation is 



p'R^l^r) 



1 d 



r'^~^ dr 

rR 

+ 2^i dr'W{r,r')u''l^^{r' 



(20) 



^ For S'-wave scattering in three spatial dimensions, ao,3 has dimensions of length and is usually called 
scattering length. 
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We rescale the radial wave function R^l\{r) as 



ut{r) = {prf^'^"RtAr\ (21) 

and obtain 



„2 (P) ( . (p) , . (2L + rf-l)(2L + rf-3) (p) 



+ 2/i / dr'W{r,r')u^ll{r'). (22) 
Jo 



The normahzation of u"l\{r) is chosen so that for r > R, 



^^^^ ^^^^(p) ^ SL+d/2-3/2ipr) + CL+d/2-3/2ipr)] • (23) 



Here Ja and are the Bessel functions of the first and second kind, Sa and Ca are the 
Riccati-Bessel functions of the first and second kind, and SL,d{p) is the phase shift for partial 
wave L. Our conventions for the Bessel functions and Riccati-Bessel functions are given in 
Appendix A In the following, we use the notation uf\{r) as shorthand for the limit p — )■ 0, 

4°)(r) = limna(r). (24) 
For our first example, d = 3 and L = 0, the rescaled radial wave function for r > R is 

Hsir) = [cot5o,3(p)^i/2(pr) - Yi/2{pr)] 

sin [pr + (5o,3(p)] 



sin [5o,3 

This satisfies the radial equation. 



(25) 



dr 

For d = 3 and L = 1, we find 



Airl(-) = -;£<l(0- (26) 



Ui,i{r) =Py^ [cot5i,3(p)J3/2(pr) - l3/2(pr)] 
sin [pr + 5i,3(p)] — pr cos \pr + 5i 3 



rsin [SiM] 

In this case 

r,^,M(r\ = ( - 



(27) 



AS(r) =(-:£ + I 1 u^Ur). (28) 
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We choose two values for the momenta pA and ps- We use the shorthand Sa = SL,d{pA), 
= ^L,d{PB)- Similarly ^^(r) = u^l^{r) and UB{r) = u^l'^\r). Therefore 

p\uA{r) = -u"Ar) + ('^L + d-mL + d-?>) ^^^^^ ^ £ r')uA{r'), (29) 

Plusir) = -u^r) + (^L + d-in2L + d-3) ^^^^^ ^ ^ dr'W{r,r')uB{r'). (30) 



We now multiply Eq. rt29| by ub-, multiply Eq. (30) by ma, and subtract the two 



{pa-Pb) UA{r)uB{r) 

= -UB{r)u\{r) + UA{r)u"ij{r) 

+ 2fi dr' [uB{r)W{r,r')uAir') — UAir)W{r,r')uB{r')] . 
Jo 

Integrating from radius p to some radius r > R, we get 

/•r 

[Pa-Pb] / dr' UA{r')uB{r') 

pr pr 

= — dr' UB{r')uA{r') + / dr' UA{r')uB{r') 
J p J p 

pR pR 

+ 2/i / dr dr' [uB{r)W{r,r')uA{r') - UA{r)W{r,r')uB{r')] 
Jp Jo 

and therefore 



(31) 



(32) 



-Pa) / dr' UA{r')uB{r') 
Jp 

= {ubu'a - uau'bWp 

pR pR 

— 2p dr dr' [uB{r)W{r,r')uA{r') — UA{r)W{r,r')uB{r')] 
Jp Jo 



(33) 



By assumption the interaction is sufficiently well-behaved at the origin and admits a 
regular solution. In particular, 



lim Ub (p) u'a (p) = lini ua (p) u'b (p) = 0. 

p-s>0+ p-s>0+ 



(34) 



So we have 



Ub (r) u'a (r) - ua (r) u'b (r) = {p% - Pa) / dr' UA{r')uB{r'). 



(35) 



The left-hand side of Eq. (35) corresponds with the Wronskian of and ua, W[ub, UA]{r), 
and so we have 

W[uB,UA\{r) = {pl-pl) [ dr'uA{r')uB{r'). (36) 



In the following, it will be useful to rearrange Eq. (23) and express u"1\{t) for r > R'm 
terms of functions r) and c(p, r) so that 



cot 6L,d{p) - 5(dmod2),0- In {PPL,d 



s{p, r) + c(p, r). 



(37) 



Later in our discussion we derive r) and c(p, r) for each of the possible cases for 2L + d, 
and show that both functions are analytic in p^. The first two series coefficients will be 
useful, 

r) = So(r) + S2(ry + O(p^), (38) 
c(p, r) = Coir) + C2(r)p2 + O(p^). (39) 
Combining with the effective range expansion, we find that 



-Soir) + Coir) + 



\rL^dSo{r) ^S2(r) + C2(r) 



p' + Oip' 



(40) 



The Wronskian W[uB,UA]{r) for r > i? is then 



W[uB,UA]ir) = (pI-pI) { ^rL4W[so,Co]ir) + 



W[s2,So]ir) 



-W[s2,Co]ir) 
+ 0{p\) + Oip%) 



-W[c2,So]ir) + W[c2,Co]{r) 



(41) 



Starting from the Wronskian integral formula Eq. (|36|), we set pa = and take the limit 

(42) 



P = P_B — ^ 0. When combined with the expansion in Eq. (41) we find that for r > R 

2 



-rL,dW[so,Co]{r) = bLAr)-'2 dr' uf\{r') 



where 



hLAr) = 2W[c2,co]{r) - — {W[s2,co]{r) + W[c2,so]{r)} + ^W[s2, so]{r). (43) 



''LA 



A fundamental result for second-order differential equations known as Abel's differential 
equation identity implies that the Wronskian of s{p,r) and c{p,r) is independent of r. 
Given our choice of normalization, the Wronskian of s{p, r) and c{p, r) is — 1 for all p. This 
implies that 

W[so,co]ir) = ~l, (44) 

W[s2,co]{r) = W[c2,so]{r). (45) 

In Appendix |A| we check explicitly that these identities hold in all cases and derive explicit 
expressions for hiAr). 

For our first example, d = ?> and L = 0, the functions s(p, r) and c(p, r) are 



s{p,r) 



smpr 
P 



(46) 
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c{p, r) = cos pr. 

The low-momentum expansions of these functions are 



s{p, r) = r- —p^ + 0{p^), 
o 



c{p,r) = 1 - — / + 0(p^). 



From Eq. (43) the Wronksians of the expansion coefficients give 

2r2 



60,3 (^) = 2r 



+ 



2r^ 



00,3 Sq-o.s 



(47) 

(48) 
(49) 

(50) 



For our second example, d = 3 and L = 1, the functions s{p,r) and c{p,r) are 

s{p,r) 



1 / sinpr 

— cospr 

p2 y py. 



. cospr . 

dp, r)=p\ h smpr 

pr 



In this case the low-momentum expansions are 

s{p,r) 



2 4 



X V 

c{p,r) = - + + O(p^). 
r 2 

The Wronksians of the coefficients lead to the result 

, , , 2 2r2 2r5 
&i,3 r = + 



r 3ai : 



45a2 ■ 



(51) 
(52) 

(53) 
(54) 

(55) 



IV. CAUSALITY BOUNDS 

We are now in the position to write down the causality bound for the effective range r^^, 



Using the Wronskian identities, Eqs. (44, 45), we can simplify Eqs. (42, 43) to 

2 

dr' u'-^^air') , 
Jo L 

hiAr) = 2W^[c2,Co](r) - —W[s2,c^]{r) + ^W[s2, So]{r). 



rL,d = bL,d{r) -2 dr' uf^\{r') 



(56) 



(57) 



"LA 



These equations hold for any r > R. In Appendix [A], we derive explicit expressions for 
the quantity &L,d(r) for all relevant combinations of d and L. In particular, we find for 
2L + ci = 2: 



2r' 



TT 



In 



1 vr 
+ 7 \ 

2 2aL, 



-\ 2 



+ 



(58) 
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for 2L + d = A: 



4 

TT 



In 



2p 



LA 



+ 7 



+ 



(59) 



and when 2L + is any positive odd integer or any even integer > 6 
hLAr) = 



2T{L + f - 2)r(L + f - 1) X -2L-d+4 

2. 



vr 

1 /rx2 



+ 



laL,d V2 
27r 1 /r\2L+'i 



T{L+l)T{L+l + l)al, 



(60) 



Since the integrand on the right-hand side of Eq. (56) is positive semi-definite, r^^d satisfies 
the upper bound 

rL4 < bLAr) (61) 



for any r > R. Eq. (61) together with Eqs. (58,|59, 60) constitutes the generahzation of the 
causahty bound on the effective range for arbitrary dimension d and angular momentum L 

We can extend our causahty bound to nonzero values of p. Consider any p ^ 0. If 
sihSlAp) 7^ 0; then W[uB,UA]{r) ioi r > R is an analytic function of pa and pb in a 
neighborhood of p. In this neighborhood, we can also consider [m^, 'u^](r) as an analytic 
function of the variables — p\ and + p\. Since W[ub, UA]{r) is antisymmetric with 
respect to pa and ps, it is an odd function of p| — Pa- Hence 



W[uB,UA]{r) 

pI-Pa 



(62) 



is also analytic with respect to p^ — p\ and + Pa- Taking the limit pb Pa, we find 
that 



W[uB,UA\ir) 



lim 

PB^PA p% - p\ 



dr' UA{r'). 



(63) 



Since the right-hand side is non-negative, we conclude that for any Pa ^ ^ such that 
sin SlApa) 7^ 0, 



lim m^p^ > 0. 

PB^PA p% - p\ 



(64) 



This is the generalization of the causality bound for nonzero p. The equivalence of our 
causality bound with Wigner's original bound on the energy derivative of the phase shift [1] 
is demonstrated in Appendix |B] 



V. CAUSALITY CONSTRAINTS ON LOW-ENERGY UNIVERSALITY 



We now discuss the impact of the causality constraints from Eqs. (61, 58, 59, 60) on 
low-energy universality. 
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FIG. 2: Hierarchy of terms in the effective range expansion as p — )• 0. 

We consider the scattering amplitude in the low-energy limit p — )• while keeping the 
interaction range R fixed. Let aL,d{p) describe an effective scattering parameter, 

The scattering amplitude is proportional to aL4{p) times a factor of p^^ from the angular 
momentum projection, 

IlAp) ^P^^o^lAp)- (66) 
In the limit p — ?■ 0, the hierarchy of terms in the effective range expansion depends on the 
value of 2L + d. This hierarchy is sketched in Fig. [2] In particular, the effective range 
parameter is as important at low energies as the unitarity contribution for 2L + = 4 and 
becomes more important for 2L + d > 5. This implies that the scale-invariant unitarity limit 
can not be reached in those cases because the Wigner bound prevents the effective range 
from being tuned to zero. In the following, we discuss the various cases in detail. 
For 2L + d = 1, we find that 

= -ip + + 0{p'). (67) 

As p — )■ 0, the effective scattering parameter has a scale-invariant weak-coupling limit 

(^lAp) ~ -iP 0' (68) 

with the dimensionful parameter, a^^^, determining the leading correction to scale-invariant 
physics. 

For 2L + d = 2,we find 

aLAp) = ^ , + OipVln'p), (69) 

2[n {-ippL^dj 
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where pL,d denotes the special value for p^^d that makes the inverse scattering parameter 



'^/0'L4 on the right-hand side of Eq. (15) equal to zero. As p — )■ 0, the effective scattering 



parameter has a logarithmic weak-coupling limit 

parameterized by the length parameter Pl4- 

For 2L-|-(i > 3, it is more convenient to consider the inverse effective scattering parameter. 
For 2L + = 3, 

«zi(p) = «zi + ^p + o(p'). (71) 

If we fine-tune the interaction so that = 0, then as p — ?■ the effective scattering 
parameter has a scale-invariant strong-coupling limit 

«Zi(p) ^^P^0■ (72) 

For d = 3 and L = 0, this is the physics of the unitarity limit in three dimensions. 
For 2L + d = 4, we have 

2 

«Li(p) = "Zi - -P"^ In {-WPL,d) + C'(p^), (73) 

where p^^d denotes the special value for p^ d that makes the dimensionless effective range 
parameter r^^d on the right-hand side of Eq. (fisj) equal to zero. After fine-tuning = 0, 



then as p — )■ the effective scattering parameter has a logarithmic strong-coupling limit 

2 

"Z ~ In {-WPL,d) ^ 0. (74) 

TT 

We note the emergence of a second relevant dimensionful parameter, pi^d- In the limit 
\(^L,d\ — ^ oo, our causality bound for 2L + = 4 places an upper bound on pL,d, 

PL,d < (75) 

For 2L + d > 5, we have 

( 0{p^) for 2L + d = 5. 
«Zi(p) = ^L4 - orL,dP' + { 0{p' \np) for 2L + d = 6. (76) 

0{p^) otherwise. 

Again we first fine-tune a~^\ = 0. This produces a power-law strong-coupling limit propor- 
tional to the dimensionful parameter, r^, ^, 

"Zi(p)~-^^L,<iP'^0- (77) 

Since this is not scale invariant we might consider a second fine-tuning where ri.d is also 
tuned to zero. However, this is not allowed by causality. In the limit \aL^d\ oo our 
causality bound for 2L + d > 5 places an upper bound on r^^d, 

< 2r(L + f-2)r(L + f-i) rny^'^' ^^^^ 



TT V 2 
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Since the expression on the right-hand side is a fixed negative number, the value r^^d = is 
not allowed. 

We see that for 2L + d>Awe are left with two relevant parameters which parametrize the 
strong-coupling low-energy limit. This corresponds to two relevant directions near a fixed 
point of the renormalization group, and the universal behavior is characterized by two low- 
energy parameters. For the case of P-wave neutron-alpha scattering in three dimensions, 
this issue was already discussed in [21]. Proper renormalization of an effective field theory 
for P-wave scattering requires the inclusion of field operators for the scattering volume and 
the effective range at leading order. In the renormalization group study of [9J , the emergence 
of two relevant directions around a fixed point was observed for various model potentials. 



VI. IMPACT OF CAUSALITY BOUNDS 



In the following, we illustrate the impact of the causality bounds for three examples. We 
start with a spherical step potential in d dimensions. This corresponds to a purely short- 
range interaction and our causality bounds apply strictly. As a second example, we consider 
the neutron-alpha interaction which is characterized by resonant P-wave interactions. The 
interaction is mediated by pion exchange which corresponds to an exponentially-bounded 
interaction of 0(e~'"/^) at large distances. In this case, the results should still be accurate 
with only exponentially small corrections. Finally, we consider the long-range van der Waals 
interaction where our general bounds do not apply. We show how our treatment must be 
modified in this case relevant to ultracold atoms. 



A. Spherical step potential in d dimensions 



As an example of the results discussed we consider a spherical step potential with radius 
R and depth V^tep, 

Wir,r') = V,,,^e{R-r)5{r-r'). (79) 

We define k so that = — 2yuV^tep and p' = \/j)^~-\-~k?. A repulsive step corresponds with 
positive imaginary k, and an attractive step corresponds with a positive real k. For the 
exterior region, r > R, the wave function u^l^^{r) is given by Eq. (23). For the interior 
region, r < R, the wave function is 



u 



(p) 

LA 



, .p'r-K 
r) (X \ l —^J, 



L+d/2-l 



{p'r) 



(80) 



Matching at the boundary r = R we find 



cot 6L,d{p) 

^ pRJajp'R) [Y^-ijpR) - Y^+i{pR)]-p'RY^ipR) [Jg^ijp'R) - Jg+ijp'R)] 

pRJaip'R) [Ja-l{pR) - Ja+l{pR)] -p'RJa{pR) [Ja-l{p'R) ' Ja+l{p' R)] ' 

where a = L + d/2 — 1. 



(81) 



We use Eq. (|81|) to generate the effective range expansion for 

2 



p 



2L+d-2 



cot 6l4{p) - (5(dmod2),o- In {ppL,d) 

n 



(82) 



14 



For 2L + d = 1, the lowest two coefficients in the effective range expansion are 

I kR + cot (kR) 



kR 



.3 _ kR {2k^R^ - 3) + 3(2K2i?2 _ 1) cot (kR) + SnRcot^ (kR) 



(83) 

In each case we multiply by powers of R to render the quantity dimensionless. For 2L+d = 2, 
we use the convention pL,d = ^- In this case the dimensionless coefficients are 



''LA 



2 

TT 



7 + 



Mf^R) 

kRJi{kR)_ 



rrdR~ 



■KH?R'^ 



\,'R'-2 + 2^R^'^''^^ 



'Ji{kR) 



For 2L + d = 3, 



_^ uRcos (kR) 

^^''^ " kR cos {kR) - sin (kR)' 



(85) 
(86) 

(87) 



rL,dR 



-1 



2k^R^ + 2kR {k^R^ - 3) cos {2kR) + 3{-2k^R^ + 1) sin {2kR) 
6k,R [— Ki?cos (kR) + sin (kR)] 

R 



For 2L + d = 4, we again use the convention pL,d — 2 ' 



4 Jo (kR) 

71 J2 (kR) 



(89) 



rL,d = {8 (27 - 1) kRJ^ (kR) + kRJI {kR) [(47 - 3) k^R^ - 4] + AkRJq (kR) J2 (kR) 
+8 Jo {i^R) Ji (i^R) [(-27 + 1) k'^R^ + 1] } X ^ 



For 2L + d = 5, 



ttkR [kRJo (kR) - 2 Ji {KR)y 
3k^R^ 



-1 d3 ^ 

^''^^ -3 + /€2i?2 + 3^^cot(fi:i?)' 



(90) 



(91) 



r^^R = {-18kR [k^R^ + 5) + ISkR [k^R^ - lO) cos {2kR) 

kR 

+45 {-2k'^R^ + 3) sin (2/?/?)} x 



9- (92) 

10 [3Ki? cos {kR) + {k^R^ - 3) sin {KR)f 

In Figs. [3] and |4| we plot the function [6i,,d(r) — 7^2L+a!-4 function of r/R for the 
sample values 2L + ci = 2, 3. We see that the function is non-negative for r/R > 1, as 
required by causality. As k^R"^ — t- —00, the potential becomes a hard spherical barrier with 
the wave function u'l^ir) vanishing in the interior region, r < R. Since 



dr' 



1 2 



^0^ 



(93) 



the causality bound is saturated in the hard barrier limit for r = R, 

biAR) - rL,d ^ 0+. (94) 
For the other values of 2L + J not shown in Figs. [3] and |4| the behavior is qualitatively the 



same. 
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rlR 

FIG. 3: Plot of \bL^d{r) — r^^^d] / B? as a function of r/R for spherical step potentials for 2L + d = 2. 




0.5 1 1.5 2 
r/R 



FIG. 4: Plot of [bL,d{r) — ri^d] /R as a function of r/R for spherical step potentials for 2L + d = 3. 
B. Alpha-neutron scattering 

Our results are exact only for the case where the interaction vanishes for r > R. For 
exponentially-bounded interactions of 0(e~''/^) at large distances, the results should still 
be accurate with only exponentially small corrections. For an exponentially-bounded but 
otherwise unknown interaction, the non-negativity condition for bL,d{r) — ri^d can be used 
to determine the minimum value for R consistent with causality. One example of an 
exponentially-bounded interaction is the three-dimensional scattering of an alpha particle 
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1 2 3 4 5 
r [fm] 



FIG. 5: Plot of bL^3{r) — r^^s as a function of r for alpha- neutron scattering in the 5*1/2 ; Pi/2: &iid 



P3/2 channels. 



and neutron. We consider S-wave and P-wave alpha-neutron scattering. 
For S'-wave scattering, 

2^2 



6o,3(r) = 2r + (95) 

0-0,3 -J^CS 



and for P-wave scattering, 



r 3ai,3 45ai3 

In Fig. [5] we plot &L,3(r) — ri,,3 for the S'1/2, P1/2, and P3/2 channels. A qualitatively similar 
plot was introduced for nucleon-nucleon scattering in the S'-wave spin-singlet channel |32] . 
We use the values 09,3 = 2.464(4) fm and ro,3 = 1.39(4) fm for ^1/2; 01,3 = -13.82(7) fm^ 
and ri,3 = -0.42(2) im^ for P1/2; and 01,3 = -62.951(3) fm^ and ri,3 = -0.882(1) fm^ for 
P3/2 [33J. The non-negativity condition for &L,3(r) — ri,3 gives R > 1.1 fm for 5*1/2, R > 2.6 
fm for Pi/2, and R > 2.1 fm for P3/2. For comparison, the alpha root-mean-square radius 
and pion Compton wavelength are both about 1.5 fm. Since the minimum values for R 
are not small when compared with these, some caution is required when choosing the cutoff 
scale for an effective theory of alpha-neutron interactions. 

We can use the results for the spherical step potential to reproduce the scattering pa- 
rameter and effective range parameter for alpha-neutron scattering. For the S'-wave, 
ro,3ao3 = 0.564. This gives k^R^ = —5.56. In turn this implies tq^s/R = 0.329 and 
therefore R = 4.23 fm. For the Pi/2-channel, r^^a^l = 0.977. This gives k'^R^ = 5.034, 
ri 3P = —1.67, and so R = 3.97 fm. For the P3/2-channel, r^^a^^l = 0.0232. This gives 
«;2p2 ^ 8 n 3P = -2.84, and R = 3.22 fm. 
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C. Van der Waals interaction 



The physics of long-range interactions must be treated separately since each long-range 
behavior determines its own low-energy universality class. For cold alkali atoms our analysis 
must be modified to take into account long-range van der Waals interactions of the type 



W{r, r') = -Cer-'^dir - r') 



(97) 

for r, r' > R. It is convenient to reexpress Cq in terms of the length scale (3q = {2fiCQY^^. 
In the following, we set c? = 3 and drop the d subscript. Instead of free Bessel functions, 
scattering states should be compared with exact solutions of the attractive potential 
[Ml [35] . The effect of the interactions for r < R are described by a finite-range i^-matrix 
Kl{p'^) that is analytic in [36] . 



(98) 



n=0,l, 



When phase shifts are measured relative to free spherical Bessel functions, the effective 
range expansion is no longer analytic in p"^. For L = 0, the leading non-analytic term is 
proportional to p^, 



pcot 6o{p) 



[r(i/4)]^/36(^r 



+ (3i + SK 



(2) 



[r(l/4)]' /^i (^( 



(0) 



GOtt 



A-r - 1 



^p^ + 0{p'^ \np). 



A-r - 1 



-p 



(99) 



For L = 1 the leading non-analytic term is proportional to p^ , 

p^ cot 6i (p) 



i8[r(3/4)]^irr mnmT[K[ 



'(0) 



71(31 



(0) 



357r/3| 



Kf^ + 1 



-p 



-4410 [r(3/4)]^ 


) + l" 


] 


/3| 






>+ 5832 [r(3/4)]^/3|(i^f))' 


12257r/3| 


'k? + 1 


3 



p 



+ 0(p^). 



(100) 



This term voids the standard definition of the effective range parameter for P-waves. 

However, one can still obtain useful information. The zero-energy resonance limit is 
reached by tuning the lowest-order i^-matrix coefficient K^l^"* to zero. In this limit the 
leading non-analytic terms in the effective range expansion vanishes, and one can define an 
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effective range parameter for both 5"- and P- waves [351 EZ] , 



[r(i/4)]^ 



To 



(51 + Sir, 



(2) 



(101) 



-36[r(3/4)]^ Pl-hK\ 



.(2) 



57TI3I 



(102) 



For the case of single-channel scattering of alkali atoms, the coefficients are negligible 
compared with /3g . This is also true for some multi-channel Feshbach resonance systems |3H] • 
In these cases we observe that the upper bounds for ri, in Eq. (60) are satisfied for L = 
and L = 1 when we naively take R ^ Pq. In general, there may be multi-channel systems 



(2) 

where the coefficients cannot be neglected. Nevertheless, 
satisfy causality bounds similar to those derived here for the effective range parameter. 



the coefficients K^j^'' should 



VII. SUMMARY AND OUTLOOK 

In this paper, we have addressed the question of universality and the constraints of 
causality on quantum scattering processes for arbitrary dimension d and arbitrary angular 
momentum L. We have derived the Wronskian identity for two solutions of the radial 
Schrodinger equation with different energy and used this identity to generalize the causality 
bound on the effective range to arbitrary d and L. Moreover, we have derived a general 
causality bound for energies away from threshold. 

For finite-range interactions, we have shown that causal wave propagation can have sig- 
nificant consequences for low-energy universality and scale invariance. For 2L + d > 4, two 
relevant low-energy parameters are required in the strong-coupling low-energy limit. In the 
language of the renormalization group, this corresponds to two relevant directions in the 
vicinity of a fixed point. In particular, we confirm earlier findings for three-dimensional 
P-wave scattering [21] based on renormalization arguments and for higher partial waves in 
general [S] in the framework of the renormalization group. 

In the low-energy limit, the hierarchy of terms in the effective range expansion depends 
on the value oi2L + d (cf. Fig. |2]). In particular, the effective range parameter is as important 
at low energies as the unitarity contribution for 2L + d = 4 and becomes more important 
for 2L + d > 5. Our results imply that the scale-invariant unitarity limit can not be reached 
in this case because the causality bound prevents the effective range from being tuned to 
zero. This has important consequences for the universal properties of systems with P-wave 
and higher partial wave interactions. Causality also constrains the wave functions and the 
probability to find particles at large spatial separation [2]. This issue will be discussed in a 
separate publication. 

We stress that our results strictly apply only to energy-independent interactions and 
single-channel systems. For energy-dependent interactions it is possible to generate any 
energy dependence for the elastic phase shifts even when the interaction W{r, r'; E) vanishes 
beyond some finite radius R for all E. Under these more general conditions, there are no 
longer any Wigner bounds and the constraints of causality seem to disappear. However, 
it is misleading to regard interactions of this more general type as having finite range. As 
noted in the introduction, the scattering time delay is given by the energy derivative of the 



19 



phase shift. The energy dependence of the interaction can by itself generate large negative 
time delays and thereby reproduce the scattering of long-range interactions. In this sense 
the range of the interaction as observed in scattering is set by the dependence of W{r,r'; E) 
on the radial coordinates r, r' as well as the energy E. In this case our bound can be 
viewed as an estimate for the minimum value of this interaction range. For coupled-channel 
dynamics without partial wave mixing the analysis can proceed by first integrating out 
higher-energy contributions to produce a single-channel effective interaction. In order to 
satisfy our condition of energy-independent interactions, this should be carried out using a 
technique such as the method of unitary transformation described in Ref. [391 - HT] . 

Our analysis concerns only the question of universality in two-body scattering. Uni- 
versality for higher few-body systems requires a detailed analysis for each system under 
consideration. For resonant S-wave interactions, the question of universality has already 
been explored for three and more particles [S]. In two dimensions, the properties of A^- 
boson droplets are universal for large, but below some critical value |12] - H^ . In three 
dimensions, the Efimov effect generates a universal spectrum of shallow three-body bound 
states |45j and two universal four-body states are attached to each three-body Efimov state 
[i6t HZj. Effective field theory and renormalization group methods may provide a useful 
starting point to extend these studies to more particles and higher angular momentum [18]- 
|5T] . Our results help to clarify some of the conceptual and calculational issues relevant to 
few-body systems for general dimension and angular momentum and their simulation using 
short-range interactions. 
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Appendix A: Expressions for 6/,d(r) 



In this Appendix, we check explicitly that the identities Eqs. (44, 45) hold for all combi 



nations of d and L and derive explicit expressions for hi^cL{r). Moreover, our conventions for 



the Riccati-Bessel and Bessel functions are given in Eqs. (A4, A5, A14, A15). 



1. Positive odd integer 2L + d 



For positive odd integer 2L + d li is convenient to use the second line of Eq. (23). For 

r>R, 

(P) (rA - „i+d/2-3/2 



pL+d/2-3/2 ^^^^ 6LAP)SL+d/2-3/2{pr) + (pr)] . (Al) 



Writing this in the form dictated in Eq. (37), we have 

sip, r) = p-^-'^/2+i/2^^+,/2_3/2(pr), (A2) 
c(p, r) = p'^'''/'-'/'CL+d/2-3/2{pr). (A3) 



20 



When 2L + d is a positive odd integer, L + (i/2 — 3/2 is an integer greater than or equal to 
— 1. For integer n, the Riccati-Bessel functions are given by 

°° •2mo — 2m— n— 1 

SJx) = v^s"+i y — — ^x^"" (A4) 

1 11°'^ „'2mr)— 2m+n, 

Cn{x) = -^x-r{-n + -)r(n + fT-Trnv Ttt^""- (^5) 

VTT / 2 L (m + 1)1 (m — n + i) 

* m=0 ^ ' ^ ^' 

The first two series coefficients for s{p,r) and c{p,r) in powers of are 



^"W-f^y • = -F(ZTiTT) (2) • 



/V V2/ , .V / ^ ^2 

The corresponding Wronskians are 

M/[so,co](r) = -l, (A8) 

TT /7-\2L+(i 

^l----°lW = r(,^|)r(L + i + i) (2) • (^^) 

2 

W[,2,co](r) = H^[c2,.o](r) = 2(2^|^_2) , (AlO) 
W^[c2,Co](r) = -iAil±^^liAfl±^^n . (All) 



TT V2 

We conclude that for r > R, 

, , , 2r(L + f - 2)r(L + f - 1) /r\-2L-d+4 
&L,d(r) = I - 

TT V/ 

4 1 /r\2 



L + |-laL,d V2 



27r 1 /r\2i+rf 

+ r(L + |)r(L + 1 + 1)4; UJ • ^ ^ 



2. Positive even integer 2L + d 



For positive even integer 2L + it is convenient to use the first line of Eq. (23). For 
r> K 



,{p) 



P^pL+d/2-3/2 ^^^^ 6LAp)JL+d/2-iipr) " Fz^+./^-i (pr)] . (A13) 
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When 2L+d is a positive even integer, L+d/2 — 1 is a non-negative integer. For non-negative 
integer n we have 

(-1)™ /X\2'"+'^ 



oo 

m=0 



T{m + l)T{m + n+l) \2 



and 



TT ^ T(m + l)T(m + n + l) \2J 



m=0 



7 is the Euler-Mascheroni constant, and is the k^^ harmonic number. 



k 

^ — ^ m 



m=l 



We can write u^I\{t) in the form dictated in Eq. (37), if we let 



s(p,r) 



>L+d/2-i(^'r) In {ppL,d) JL+d/2-i{pr) 

TT 



We now consider each of the possible cases for positive even integer 2L + d. 



(A14) 



(A15) 



(A16) 



(A17) 
(A18) 



a. Case 2L + d = 2 
When 2L + (i = 2 we have 

c(p, r) = - 




rir 




Yo{pr) In {ppL,d) Mpr) 

TT 



Both functions are analytic in p^. The first two series coefficients are 



soir) 



r\i/2 



Coir) 













+ 7 







^ \2. 

The corresponding Wronskians are 



C2(r) = 



r-N 5/2 

"2) 



2 /r\5/2 



In 



2p 



W[so,co]ir) 



(A19) 
(A20) 

(A21) 

7-1 
(A22) 

(A23) 
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W[s2,So]{r) = TV 



W[s2,co]{r) ^W[c2,so]ir) ^ -2 
4 



W[c2,ccl(r) 



TT 



+ 



1 



The function bL,d{r) for this case is 

2r2 



TT 



In 



+ 7 



1 TT 



2 2aLw 



+ 



b. Case 2L + d = A 
For 2L + d = 4, 

c{p, r) 

The series coefficients are 



and the Wronskians are 



s(p,r) 



7rr 





P Ji{pr), 



Yi(pr) In {ppL,d) Jiipr) 

TT 



TT /r 



2 V2 



7/2 



C2(r) = 



1 /r\3/2 



V2 

iy[so,co](r) = -l, 



2 In 



+ 27 



TT /r\4 
2 V2/ 



W^[s2,so](r) 
W^h,co](r) = l^[c2,so](r 



r\ 2 
2. 



W^[c2,Co](r) 



TT 



In 



+ 7 



The function hL,d{j') is 



4 

TT 



In 



+ 7 



4 /r\2 TT /r\4 
+ ^ 



aL,d ^2 



^2. 
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c. Case 2L + d>6 



The last case we consider is when 2L + is an even integer greater than or equal to 6. 
Here we have 

sip,r) = ^p-^-'i/^+\j^_^,/,^,(pr) (A37) 

c{p, r) = -y^Z+'^/^-i YL^,/2-i{pr) - ^ In (pp^,,) JL+d/2-i{pr) • (A38) 

The first two series coefficients have exactly the same form as in the case for positive odd 
integer 2L + d, 



go(r)= ^\ (-) , S2(r) = ^ (-) , (A39) 

r(L+ ^ - 1) /r\-i-rf/2+3/2 r(L+ ^ -2) /r\-i'-'^/2+7/2 

^ (2) • ^-w= (2) ■ '^^"^ 

We conclude the same result for &L,d(r) as for odd 2L + c/ as written in Eq. (A12). 
Appendix B: Equivalence with Wigner's bound 

In this Appendix, we demonstrate the equivalence of our causality bound with Wigner's 
original bound on the energy derivative of the phase shift 

Let lA{f) be a free incoming radial wave for momentum pA-, 



lA{r) = \^^^^^Y~P^a'^^'^ ^'^ [-i ■ JL+d/2~l{PAr) - >L+d/2-l(PA?^)] • (Bl) 

We are using the same phase convention as Wigner but a different normalization. From 
Abel's differential equation identity, the Wronskian of Ia and is independent of r. For 
our chosen normalization of the incoming wave, 

lA{r)n{r) - rA{r)Ur) = 2^pf^'-\ (B2) 

For r > R, 

We define aA{r) as the reciprocal logarithmic derivative of ua (r), 

^ UAirl ^ lAjr) - e'^'-I*Air) 



Then 

I\{r) - aA{r)I'X{ry 



(B5) 



e^^^- [I\{t) - aA{r)I'X{r)] = Ia^t) - aA{r)Ur). (B6) 
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We place a dot on top of a function to indicate the derivative with respect to pa- Differen- 
tiating Eq. (B6) with respect to pa, we get 



2te'''^6A [Ilir) - aA{r)mr)] + e'^'-^ [/l(r) - aA{r)i'X{r) - aA{r)I'X{r) 
= iA{r)-aA{r)i'A{r)-aA{r)rA{r). 

Solving for 5a gives 



where 



and 



FaW) 



= FA{r) + GA{r)aA{r), 
I i\{r) - aA{r)i'2{r) - e'^'^-^ [/A(r) - aA{r)i'A{r) 



2i 



GA(r] 



I\{r) - aA{r)I'X{r) 

1 lAir)IXir) - PAr)IXir) 
2i \lA{r)-aA{r)I'A{r)? 



(B7) 
(B8) 

(B9) 
(BIO) 



Both Fa and Ga can be simplified further. We replace a a using Eq. (B4) and find 



FA{r) = - 



1 Re |/A(r) [rX{r) - e-''''^IX{r)\ - P^r) [l\{r) - e-^^'^l 
~i I\{r)rj,{r)-lA{r)I'X{r) 



(Bll) 



The Wronskian identity, Eq. ( B2 ) , leads to 
1 



Fat) 



2L+d-2 



Re { lA{r)I'X{r) - Wllir) - e'^^'- lA{r)lX{r) - /;i(r)/^(r) 



2L+d-2 

Pa 



\lA{r)-aA{T)IX{T)\ 



With the same Wronskian identity, Ga simplifies to 

GA{r 

We note that 

d^(r) = lim 



2 • 



(B12) 
(B13) 



UB{r) _ UA{r) 



lim 



W[uB,UA]{r) 



PB^PA pb - Pa 



PB^PA {pb - Pa) u'j^{r)uX{r) 



(B14) 



For any pa 7^ 0, we use Eq. (63) to get 



aA{r) 



2p, 



dr' [ua {r')f ■ 



WAir)f 

We see that both GA{r) and dyi(r) are non- negative and so 

SA = FA{r) + GA{r)aAir)>FAir). 



(B15) 



(B16) 



This inequality holds for all 7^ 0, and therefore also holds in the limit pa 0. This is 
Wigner's causality bound generalized to arbitrary dimension d. 
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Away from threshold, the equivalence between Wigner's bound and Eq. (64) is clear from 
the biconditional statement, 

6a > F4r) ^ a4r) = Hm , ^^"^'^f^l,, > 0. (B17) 

PB^PA {pb - Pa) u'^{r)u'^{r) 

For Pa = 0, the equivalence with our causality bounds follows from the low-energy expansion 
for W[ub,ua], 

W[ub,uaW)=pI |lr^,iy[so,Co](r) + (^J-^ W[s2,So]{r) 

- — W[s2,Co]ir) - -^W[c2,So]{r) + W[c,,Co]{r)] + 0{p%). (B18) 
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